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We consider a network of interacting resonators and analyze the physical ingredients that enable 
the emergence of relaxation-free and decoherence-free subspaces. We investigate two different situ- 
ations: i) when the whole network interacts with a common reservoir and ii) when each resonator, 
strongly coupled to each other, interacts with its own reservoir. Our main result is that both sub- 
spaces are generated when all the resonators couple with the same group of reservoir modes, thus 
\^ , building up a correlation (among these modes), which has the potential to shield particular network 

• states against relaxation and/or decoherence. 

o : 

The search for mechanisms to bypass decoherence, a subject of major concern for quantum information processing, 
has deepened our understanding of open quantum systems and led to ingenious schemes for coherence control, which 
go far beyond the quest for conditions that weaken the system-reservoir coupling 0, Q • The myriad contributions 
to this subject started, inspired by classical error-correcting codes, with quantum coding schemes for information 
\^ \ stored in a quantum memory Q. On the assumption that the decoherence process acts independently on each of 
CN ■ the qubits stored in a memory, a particular qubit encoded in a block of ancillary qubits is able to withstand a 
substantial degree of interaction with the reservoir without degradation of its information. Another strategy, the 
so-called engineering reservoirs 0, compels the system of interest, whose state is to be protected against decoherence, 
to engage in additional interactions besides that with the reservoir. This program, based on the indirect control of 
system-reservoir dynamics, has been developed for trapped ions and atomic two-level systems 0, 11] as well. 
■ Finally, the process of collective decoherence, where a composite system interacts with a common reservoir, has also 
OO instigated several interesting results related to what has been called a decoherence-free subspace (DFS) '9', 10, TT]. ft 
is noteworthy that while quantum error-correcting codes presuppose quantum systems that decohere independently, 
DPS - as it has been understood until the present study - is generated by distinct quantum systems coupled to a 
common reservoir. 

' In this contribution we are concerned with collective dissipation and decoherence in a network of N coupled 
Qh! resonators. We analyze the physical ingredients ruling the emergence of a DPS and, in particular, a relaxation- 
4L ■ free subspace (RPS) composed of states protected against both dissipation and decoherence, demonstrating that a 
DPS contains a RPS. We first analyze the situation, treated in the literature to date, where all the resonators are 
coupled to a common reservoir. However, as the scenario of a common reservoir is in practice rather unusual, we next 
analyze the situation where each resonator interacts with its own reservoir, which seems to be more appropriate for 
most physical systems. In the domain of cavity quantum electrodynamics, distinct reservoirs must be considered for 
distinguishable cavities, even if they present equal quality factors, inasmuch as there are no correlations, whatsoever, 
between the reservoirs. The same applies to distinguishable trapped ions or a traveling field reaching distinguishable 
5J] optical elements. There are a few particular situations where a set of quantum systems do interact with a common 
reservoir, such as an atomic sample or distinct fields inside a perfect cavity. In the former case, different atomic 
transitions couple with different reservoir modes (KM) and again, the correlations between these KM define either a 
common or distinct reservoirs, as will be demonstrated below. In the latter, we also show that the proximity of the 
distinct field modes sets the strength of the correlation function between the RM which governs the emergence of 
both RPS and DFS. 

A common reservoir. Consider the Hamiltonian H = Hs + -\- Hj, for the system of N interacting resonators 
(Hs), the reservoir (Hf/), and the interaction between the resonators and the reservoir {Hi). Assuming, from here 
on, that the subscripts m (m') and n (n'), labelling the resonators, run from 1 to N , the Hamiltonian is given by 
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i/fl = V UJkblbk, Hi^Y] Vrnk (arnbl + h.c.) , (lb) 

where (flm) is the creation (annihilation) operator of mode w™, while and 6^, are the analogous operators 
of the fcth bath mode, whose corresponding frequency and coupling strength to the mth resonator are uju and Knfc, 
respectively. We assume that the interactions between the resonators of frequencies uim and ujn are within the rotating 
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wave approximation with coupHng strength Xmn- In what follows we analyze both the weak and strong coupling 
regimes between the resonators, where iVAm„ ^ and NX^an ^ i^mh respectively p^ . We also investigate the 
special case of the weak coupling regime, of negligible mutual interaction between the resonators, where Xmn ^ 

{uJrn' — ^n')- 

The Hamiltonian Hs can be brought into a diagonal form Hs = ^m^ni^m through a linear transformation 
Am — X^n^™"*^"' {^m} being the normal modes of the network and Cmn the elements of an orthogonal and 
symmetric matrix (C,„„ = Cnm = Cmli)- The interaction goes into the form Hj = J2m.n,k^"inVmk (^^nbl + h.c.^ 

Therefore, instead of N coupled resonators interacting with a common reservoir, as in Eqs. lla() and Ijlbjl. the 
transformed Hamiltonian describes N decoupled resonators under the action of the reservoir. In the picture defined 
by the transformation cr''-(^s+HR)t ^ obtain the interaction V(i) ~ X)m « (^mji(^)^n +h.c.), where C'|„„(t) — 

X]fe ^nin^fc^l analogy with the developments in Refs. and noting that Cmn < 1, we obtain, at 

T = OK, the master equation of the reduced system of resonators 

'^'''-f^*^ ^l[pi,...Mt),Hs]+Y^ Cmn{t)pi,...Mt), (2) 
at — 'm.n 

where the influence of the reservoirs is represented by the Liouville operators Cmn', accounting for direct (m = n) and 
cross-decay (m ^ n) channels |12| . which link together the many resonators via a common reservoir. The Liouville 
operators are given by 

>Cm„(i)pi,...,jv(i) = { [anPi,...,N{t),alJ^ + h.c.} . (3) 

where the coefficients ^mn = X)m' n> ^mm' i^n') Cm'n'Cn'n, take into account the elements of the transformed matrix, 
coming from the topology of the network: the frequency distribution of the resonators, their coupling pattern and 
strengths, ^mn also depends on the interaction between the system and the reservoir through the function Smn i^m') = 

Jo ir ('^) ^"^ ('^) ('^) e-^("-"-'H*-*') , where the frequencies of the reservoir are assumed to be continuous 
variable with spectral density cr(w). The function {^m') is related to the correlation between the RM through 

fodt' iOmm'{t)Ol^,{t')) « \emn (17„0 e^'^""' . 

To analyze the behavior of the correlation Smn i^m') we first assume, as usual, that both the spectral density a {ui) 
and the coupling strength Vmijjj) are functions that vary slowly around the frequency fim, such that Smn {^m') ~ 
a (f^m') Vm (Om/ ) Vp . (Jim')- The systcm-rescrvoir interaction, which applies to all coupling regimes, is thus modeled 
by the function [12] 

Vm (iv) = Vm e-«^(^-"^)'^ ^ . (4) 



We observe from Eq. (0J that when the separation between two normal modes is large enough compared to the 
Gaussian widths, |rij — fij/j ^ these modes are considered distinct. The difference between the weak 

and strong coupling regime is that, in the former, the normal modes Vim are approximately equal to the natural 
frequencies uim, while in the latter, they may be quite different, depending on the coupling strength Xmn 'l^ (see 
below). In both cases the interactions between the resonators couple them with all groups of bath modes around 
all the M (< N) distinct normal modes (with j = 1, M) defined by the topology of the network In fact, 
it has been demonstrated [12| that the states of two distinct resonators, say m and n, arc interchanged many times 
before the relaxation takes place, in both weak and strong coupling regimes, the recurrence time of the states being 
approximately {N Xmn)~^ ■ This interchange of excitation is the mechanism behind the coupling of each resonator to 
all groups of RM around all normal modes. When considering resonators with negligible mutual interaction, such 
that the mth resonator practically remains at its natural frequency ujm-, expression degenerates to the expected 
result Vm (uj) = Vm e"*'"^"""") where the mth resonator couples only with the group of RM around LUm- In fact, 
for negligible mutual interaction the cross-decay channels are exceedingly small, around Xmn I iy>m' — ^n')i leading to 
recurrence times longer than the relaxation of the network. 

Under the above considerations, for both weak and strong coupling regimes, we obtain the result 
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(17,) « a (17,) VmVn e-i^i^^-^.r-iA'^^-'h') ) (5) 



revealing that Emni^m') ~ cr (f^m') l^m^l = T whenever \Vn\ « \Vm\ = \V\ and a {Qm) is a constant a, 
which is a reasonable assumption for a Markovian white noise, where T = a \V\^ . Note that only the term 
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: i contributes significantly to each sum in Eq. in both weak and strong coupling regimes. The re- 
(yim') ~ r implies that 7m„ = T and, consequently, both direct and cross-decay channels have the same 
order of magnitude, a crucial condition, as demonstrated below, for particular classes of states to be shielded against 
relaxation and/or decoherence. In the case of negligible mutual interactions between the resonators, we obtain 



(w,„') W uVmVn e 
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' , which, differently from Eq. ij^l, varies from to F, depending on 
the proximity between the natural frequencies LUm- The correlation function emn (i^m') attains the maximum value F, 
bringing about a DFS, only in the degenerate case where all the resonators have the same natural frequency. We stress 
that the maximum correlation function (emn = F) between the RM will be induced by the network itself when all its 
resonators interact with the same groups of RM — those around the degenerate frequency in the case of negligible 
mutual interactions. 

In order to characterize the RFS and DFS, from here on we consider the degenerate case where ujm — and 
Amn = A, leading to two distinct normal modes f2+ = uj + (N — 1) A and fi_ = uj — \ We also consider a 

Markovian white noise reservoir, where Smn ~ £ and Tm ~ F, and the initial entanglement between the N resonators 
is given by \^i....^n) = J^Y^t^i K/?™})- This collective state comprehends a superposition of J product states, 
each one consisting of a product of TV distinct coherent states (3^, A/" being a normalization factor . We also redefine 
the correlation function e, introducing the dimensionless parameter e = e/F varying from zero to unity. Under these 
assumptions and following the reasoning in Ref. [l^. we obtain the density operator for both regimes e = 1 and 
e< 1: 
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where 

(i_,)rt/2 



-(eNr/2+in+) 



{t)Pn ^iid the function accounting for the excitation decay is given by Qmn{t) = 

N V- ' (iV^™„-l)e-^f^-*). 

Decoherence Free Subspace. To analyze the role played by the excitation decay function Omn(^) hi the emergence 
of a DFS, let us consider a particular family of initial states |^'i....,jv) whose representative element is 



R,S 



= A4 






(7) 



where R (S) indicates the number of resonators in the coherent state a (—a) in the first term of the superposition and 
—a (a) in its second term. The remaining N — R — S resonators are in the coherent states {77}. As we are considering 
a symmetric network where all the resonators have the same damping rate, they are indistinguishable. Therefore, 
swapping the states of any two resonators m and n we obtain a state which is completely equivalent to Eq. |(7J). For 
the states defined by Eq. 10 and the evolution 10), we obtain 
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with an effective lower 



R=S=o 



with 



from which two distinct classes of states follow: the superposition states |\l/i,...,iv)/j^5^7v/2 
decay rate (compared to F) F| = (1 — e)F, since = e"'"^*/^"*^"*, and the product states |^'i,...,Ar) 
an enhanced decay rate F-f = [1 + (N — 1) e] F, since = ,y g-rTt/2-»s^+t Pqj- decoherence time of the state (0) 
we obtain 
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1*1,.. ..Jv)) = 



2|a| {R-SfT+ {R + S)-{R-S) 



(9) 



to be analyzed for both classes of states l^"!^ 



In the regime e = 1 
the condition R — S 



R.=S=N/2 



and I'i'i, 



i?,=S=0 



under both regimes e = 1 and e <C 1. 



such that F| = and F| = A^F, both classes of states belong to a larger DFS defined by 



with F| = 0, are effectively uncoupled 
R=s=o dissipate at the maximum rate 



However, while the superpositions |*i,...,Ar/^^5^^y2 
from the reservoir, defining a RFS, it is worth noting that the states |5'i,...,jv 

F| = A^F. The intermediate class of decoherence-free states, where R = S < N/2, also dissipate at the maximum rate 
F| ^ NT, leading to the conclusion that the DFS contains the RFS. The correlation e = 1 is achieved independently 
of the spectral shape of the network for both weak and strong coupling regimes; in the case of negligible interaction 
between the resonators the regime e = 1 occurs only when all the resonators have approximately the same natural 
frequency. These three situations, where the resonators interact with the same group of RM, trigger a kind of 
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phase matching, between the resonators themselves together with the RM, which is unequivocally produced by the 
cross-decay channels, as discussed below. 

For the class of states not belonging to the DFS, i.e., with R ^ S, a.s for example |^'i....,Ar)7^_i — 
N± ± |— <Xi (where a Schrodinger-cat like state is prepared in resonator 1 while all the remaining 

j resonators are prepared in the coherent states 77), we get the expected decoherence time tjj ^|^'i,....Ar)^_j^ s-o) ~ 
r~^/(2 |q;i|^), even when e = 1. Finally, in the regime e ^ 1, only the product states \'^i,...,n) ^ 5^0 remain in the 
DFS, the decoherence time for the superpositions |'I'i,...,Jv)fl^5^7v/2 increasing as the correlation e decreases. 

It is worth noting that the necessary and sufficient condition for a subspace to be decoherence- free can be established 
through the Lindblad master equation 



dpi,...,N{t) 
dt 



i[pi....,N{t),Hs] + ly] r™„([A™pi A,(t),At] +h.c.) , (10) 



in that all these subspace states {^'} must be degenerate eigenstates of all the Lindblad operators {A™}, i.e., A^ l^*) = 
c„i l^f), V m, JJJ. This is exactly what happens to Eq. (jSJ when considering, for a reservoir at T = OK, the above- 
mentioned three cases where e = 1. In fact, under this condition, Eq. (j^J can be written as 

^^ii^ = ^[p^,..Mt), Hs] + ^ ([Api,....^(i), At] + h.c.) , (11) 

with the eigenvalue equation A |\E'i^...^jv)^_g = (^r]{N — 2R) /VN^ \^i....,n) j^-g for the operator A = J2n o.n/VN- As 

far as A |5'i,...,Af)^^5^jv/2 ^ ^ while A |^'i,...,Ar)^^5^Q — rj^/N \'^i,...,n) j^-g^Q, we conclude that: while decoherence- 
free states are characterized by an eigenvalue equation, dissipative-free states follow from the fact that their associated 
eigenvalues are null 1^. To summarize, we note that the crucial ingredient for establishing a DFS is the action of 
cross-decay channels Lmn{t)» which always take place for a common reservoir. When the cross-decay channels are 
of the same order of magnitude as the direct-decay channels Lmn{t)», such that e = 1, a DFS emerges for the states 
where R — S. When R^ S = N/2, these states compose a RES. 

Strong coupling and distinct reservoirs. To get round the necessity of a common reservoir for the emergence of 
a DFS, we briefly analyze the case where each resonator is coupled to a particular reservoir. Now, however, they 
must be strongly coupled to each other for the accomplishment of the conditions leading to a DFS. Considering the 
notation employed in Eqs. (|la|l and l)lb|l . the system Hamiltonian in the degenerate case ujm ~ oJo and A„j„ = A, is 
given by 12] 

Hs^uJoV] ai,am + ^y] . (at„a„ h.c.) , (12a) 
Hu = y^^ , ujmkbl^k^mk, Hi = y^ V,nk (alnKik + li-c.) , (12b) 



where tJo = jl -I- [{N — 1) (A/2a'o)] |, under the assumption that the energy spectrum has a lower bound [l3|- As 
shown in Ref . Il2!| , in the strong coupling limit NX > cDq the master equation reads 



^^^i^#^ = *bi,...,iv(t),i?s] + V /:™„Pi,...,iv(<), (13) 

at '■ — 'ra,n 

where Hg = wo Z)™ aJnOm + (A/2) Y.ra,n{m.i.n) ("JnOn + h-c) , and the Liouville operators C.mPi,...,N{t) are given by 

AnnPl,...,Ar(0 = {(Ki/3l,...,A'(0)ani] +ll-C.)}. (14) 

The coefficient = J2n' ^mli'^rn {^n')Cn'n follows from jrnn assuming £mm' (^K') = E m {Qn') Srn77i'- Eor a sym- 
metric network, in the degenerate case described by the Hamiltonian in Eqs. (|12af) and (jl2b|l . only two distinct normal 
modes arise: fl+ = lvq + {N — 1) A and fl- = uo — X. Consequently, the damping rate E™ splits into two different 
values E+ and F~, around ^1+ and ri_, respectively, as discussed in detail in Ref. 12], where we also observe that 
7„i„ oc (E+ — E~ ) /N while oc [E+ + (A^ — 1) E^J /N. Therefore, for Markovian white noise reservoirs, where 
the density of RM are the same around both normal modes, such that E+ « E~ , the cross-decay channels disappear 
together with the possibility of a DFS for superposition states of the form l^'i,...,Ar) in both weak and strong 
coupling regime. In the weak coupling regime, we also get the cross decay channels around zero even for reservoirs 
with spectral densities distinct from the Markovian white noise, since in this case fl^ « f2_ « loq and again E+ « E~ . 
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Therefore, differently from the former case of a common reservoir, here the cross-decay channels linking the resonators 
to all reservoirs besides their own, arise only from the strong coupling between the resonators for non-Markovian white 
noise. 

When distinct reservoirs are considered, the effectiveness of these cross-decay channels depends on the condition 
NX/uq > 1 (pulling the normal modes apart from each other) instead of the correlation functions between the reservoir 
operators. Therefore, considering the strong coupling regime and appropriate spectral density of the reservoirs, such 
that ./Vr~ <C r+ and, consequently, 7mn ~ 7mm, the cross-decay channels become of the same order of magnitude 
as direct-decay channels, inducing a DFS for superpositions of the form |5'i,...,jv)/j^s- In this case we get a master 
equation similar to Eq. Hll|l . apart from the system-reservoir coupling strength NT /2, which becomes r+/2 (assuming 
that all reservoirs have the same damping rate = T). (Note the lack of the factor in the latter coupling strength, 
following from the dismissing of coupling around the normal mode $7^.) Therefore, for strongly coupled resonators, 
each one interacting with a different reservoir, we get to same decoherence-free subspace states \^i,...,n) ^-g obtained 
in the case of a common reservoir. We finally mention that the subject of distinct reservoirs have analogies with 
previous work in literature , where external control of quantum systems plays similar role as the strong coupling 
between distinct resonators. 
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